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Abstract. We present plane-parallel equilibrium models 
of molecular clumps that are supported by Alfven waves 
damped by the linear process of ion-neutral friction. We 
used a WKB approximation to treat the inward propaga- 
tion of waves and adopted a realistic ionization structure 
influenced by dissociation and ionization due to photons 
of external origin. The model clumps are larger and less 
centrally condensed than those obtained for an assumed 
ionization structure, used in some previous studies, that 
is more appropriate for dark regions. 
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1. Introduction 

Giant molecular cloud complexes (GMCs), the birth 
places of stars, are typically many tens of parsecs in lin- 
ear extent and have masses fr om 10*^ to lO^Mp) and tem- 



peratures of 10-30 K (see Hartquist et al. 1998 for a re- 
cent review). Observations of CO emission f rom GMCs 
( [Blitz fc Thaddeus 198C| ; [Wilhams et al. 1995[ ) show them 
to be composed of many smaller clumps that are a few 
parsecs in extent and contain ^ lO*^ Mq. 

The widths of CO emission lines originating in indi- 
vidual clumps are supersonic and have been attributed to 
the presence of Alfven waves having subAlfvenic velocity 
amplitudes (Arons & Max 1975). The Alfven waves con- 
tribute to the support of a clump along the direction of 
the large-scale magnetic field; the damping of the waves 
affects the degree of support that they provide. An im- 
portant and well understood mechanism for the damping 
of linear Alfven waves in a partially ionized medium is 
that due to ion-neutral friction which depends on the ion- 
ization structure ( Kulsrud fc Pearce 19691 ). Ruffle et al. 
(1998, hereafter R98) and Hartquist et al. (1993) have em- 
phasised that the dependence of the ionization structure 
on total visual extinction, should greatly influence 
the density profiles of clumps if Alfven waves contribute 
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to their support. To quantify the assertion of Hartquist 
et al. (1993) and R98, we present in this paper models 
of plane-parallel, wave-supported CMC clumps like those 
identified in the work of Williams et al. (1995), who made a 
detailed analysis of the CO maps of the Rosette Molecular 
Cloud (RMC), identifying more than 70 clumps. The mod- 
els that we have constructed are for RMC-type clumps in 
equilibrium, a restriction justified by the fact that clear 
spectral signatures of collapse have been found only when 
much smaller scale features have been resolved (see, e.g., 
Hartquist et al. 199^ ). 



We have adopted a WKB description of the wave prop- 
agation as did Martin et al. (1997) in their work on wave- 
supported clumps. Their work differs substantially from 
ours in that they used an ionization structure appropri- 
ate for dark regions. Also, we have considered inwardly 
rather than outwardly propagating waves, as many of the 
clumps mapped by Williams et al. (1995) do not contain 
detected stars and may have no internal means of gener- 
ating waves. Indeed, the waves may be produced at the 
surface of a clump by its interaction with an interclump 
medium. 

Other authors have addressed the importance of pho- 
toabsorption for the effects that the ionization structure 
will have upon a clump's dynamics. These authors have 
been concerned primarily with dense cores and/or en- 
velopes around them; cores are much smaller-scale objects 
than the clumps identified in Williams et al. (1995). Mc- 
Kee (1989) addressed the possibility that collapse in a 
system of dense cores is a self-regulating process due to 
the ionization of metals such as Magnesium and Sodium 
by photons emitted by stars formed in the collapse; he 
was concerned with infall due to ambipolar diffusion of a 
large-scale magnetic field. Ciolek & Mouschovias (1995) 
have shown that the large-scale magnetic field can sup- 
port a photoionized envelope around a dense core for a 
time that is very long compared to the ambipolar diffu- 
sion timescale in the center of the dense core. In contrast to 
McKee (1989) and Ciolek & Mouschovias (1995), Myers 
& Lazarian (1998) addressed the effect of photoabsorp- 
tion on support by waves rather than by the large-scale 
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magnetic field. They stressed that observed infall of dense 
core envelopes is slower than that expected due to gravi- 
tational free-fall and more rapid than collapse due to the 
reduction by ambipolar diffusion of support by an ordered 
large-scale magnetic field. They considered collapse of ma- 
terial supported primarily by waves and subjected to an 
external radiation field. While they made clear comments 
about the importance of the Ay dependence of the ioniza- 
tion structure for their model, they did not perform any 
detailed calculations in which a realistic dependence of the 
ionization fraction on Ay was used. 

Several sets of authors have considered nonlinear ef- 
fects in the dissipation of waves supporting a clump. 
Gammie & Ostriker (1996) investigated models of plane- 
parallel clumps and from their "1 2/3-dimcnsional" mod- 
els found dissipation times due to nonlinear effects to be 
longer than the Alfven crossing times for a fairly large 
range of parameters. The three dimensional investigations 
of Mac Low et al. (1998) and Stone et al. (1998) suggest 
the more restrictive condition that the angular frequency 
of the longest waves be no more than a few times 27r/tA 
(where is the Alfven crossing time) in order for the dis- 
sipation timescale due to nonlinear damping to be roughly 
the Alfven crossing time or more. In most cases addressed 
in this paper we have restricted our attention to such an- 
gular frequencies so that we are justified to lowest order 
in focusing on only the damping due to ion-neutral fric- 
tion. It should be noted that the above three dimensional 
studies of nonlinear effects concerned homogeneous tur- 
bulence and did not include ion-neutral damping for a 
realistic ionization structure. If we are correct in suppos- 
ing that the waves in clumps are driven externally, then 
the turbulence is not homogeneous and its nature depends 
on both the viscous scale set by ion-neutral damping and 
the exact boundary conditions. The effects of nonlinear 
damping and multiple dimensions will be considered in 
subsequent work. 

In Sect. 2 we present the equations for the wave en- 
ergy, the static equilibrium clump structure, and the grav- 
itational field. In Sect. 3 we give a description of the cal- 
culations of the ionization structure for various values of 
the clump density and Ay while Sect. 4 contains details of 
the models considered here. Finally, in Sect. 5, we present 
conclusions. 



2. Equations of Wave Propagation and Static 
Equilibrium 

We consider plane-parallel clumps with z = correspond- 
ing to the clump midplane and z = (with Zb defined as 
positive) corresponding to boundary between the clump 
and the interclump medium. The large-scale magnetic 
field is taken to be Bqz with z normal to the surface of a 
plane-parallel clump. We study waves of angular frequency 
Lo propagating from z = Zb in the —z direction. 



The ion velocity can be expressed as 
Vi = Ve J , 



(1) 



where V is defined below and is the real component 
of the complex wave vector. Hartquist & MorfiU (1984) 
used a two-fluid treatment to examine a related problem 
and showed that for an inwardly propagating linear Alfven 
wave 



d Vo 

az v^. 



(2) 



where v* is the complex conjugate of the ion velocity per- 
turbation and 



v. 



(3) 



where pi is ion mass density and is the neutral mass 
density. The ion Alfven velocity is given by 

Bo 



(4) 



The ion-neutral coupling frequency, voPii is such that the 
momentum transfer per unit volume per unit time from 
ions to neutrals is given by i^oPriPiivi — Vn) where Un is the 
velocity of the neutrals. To a reasonably good approxima- 
tion, is the dominant ion and we may take 



1^0 Pi 



2.1 X 10~^sec 



Hi 



(5) 



where rii is the ion number density ( Osterbrock 1961 ). If 
the dominant ion species are very massive, as occurs at 
large Ay and densities, the constant would approach 2.3 x 
10~^s~^. However, in this work we ignore this dependence 
and use 12 mjj as the mass per ion. 

If V in Eq. [l] is written as e J ' ^ , where fcj is the imag- 
inary component of the wave vector, and the linearized 
version of the induction equation (cf. Eq. 3 of Hartquist 
& MorfiU 1984) is used, Eq. | yields 



d 

dz \ fc? 



k 



■bb* = 



V2 bb* 
vl, fc? + kf 



(6) 



where b is the perturbation magnetic field and b* is its 
complex conjugate. In the WKB approximation, (PV/dz'^ 
is taken to be equal to zero, which is equivalent to assum- 
ing 

kf + dfci/dz < fc2 (7) 

Then it follows (cf. Eq. 7b of Hartquist & MorfiU 1984) 
that in the WKB approximation 



J^ + Vi 



vith 



Vl 



I'QUj'^PnPi 

vlpf + uj^pI 



(8) 



(9) 



Coker, Rae, & Hartquist; Turbulent Wave Support of Clumps 



3 



As is consistent with the WKB approximation, we assume 
kf <^ and substitute Eq. ^ into Eq. ^ to find 

d_ f VA,U \ V2 
dz V V^^^ + 



Table 1. Assumed CO Photodissociation Rates 



-u , 



(10) 



ViJ uj'2 + Vi 

where U — bb* /16tt is the time-averaged energy density of 
the perturbation magnetic field. 

We solve Eq. |l^ along with the static equilibrium equa- 
tion 
^d{pn + Pi) 



dz 



dU 

^ = -(Pn+Pi).9 



( Martin et al. 1997 ) and the gravitational equation 



dz 



47rG(p„ 



(11) 



(12) 



where Cg, g, and G are the isothermal sound speed, the 
strength of the gravitational field, and the gravitational 
constant, respectively. 

We verify the assumption that kf ^ kf a posteri by 
checking that 

,2 ■ 



1 dV^ 



(13) 



Note that if Eq. |l^ is satisfied, Eq. ^ is as well. 

3. Calculations of the Ionization Structure 

The ionization structure determined by R98 and presented 
in their Fig. 1 was calculated on the assumption that, due 
to shielding of the CO by itself and by H2, the rate of CO 
dissociation by photons of external origin is negligible. For 
a plane-parallel semi-infinite cloud with constant Hydro- 
gen nucleus number density, nn = 10^ cm~^, we assume 
an Ay-dependent dissociation rate that results in a CO 
abundance relative to rifj, x(CO), that is in harmony with 
the measurements shown in Fig. 6 of van Dishoeck (1998). 
For the nn = 10"^ cm"'^ model. Table [l] gives x(CO) and 
the photodissociation rate as a function of Ay. Note that 
the total abundance of carbon nuclei relative to rifj is fixed 
at 10"^. In the work reported here, we used an ionization 
fraction that depends on both Ay and nn. Using the Ay- 
dependent CO photodissociation rate from Table ^ wc 
calculated the fractional ionization as a function of Ay 
for riH = 3 X 10^,3 x 10^, and 10^ cm~'^ at various Ay 
values. A bilinear interpolation in Ay and nu is used to 
find the actual ionization fraction used for a given point in 
the clump. For densities above 10* cm"'^ we assume that 
the total ionization fraction, ^ = Pi/pn, goes as ny^^"^ , as 
expected in a dark region. Note that in all models. Ay is 
always greater than 4 whenever nn > 10''' cm~'^. 

Since the depletions in RMC-type clumps are very un- 
certain, we present results for both depletion cases given 
in R98. As discussed therein (also see Shalabiea & Green- 
berg 1995), case A abundances resemble those seen in dark 
cores with Ay ^ 5 while case B, with higher fractional 
abundances of lower ionization potential elements, is more 
appropriate for more diffuse clouds. 
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4. Details of the Models 

Many RMC-type clumps are not bound by their own 
gravity and must be confined by interclump media 



(Bertoldi & McKce 1992). We shall assume that an in- 
terclump medium is sufficiently tenuous that it does not 
shield the clump from the standard interstellar back- 
ground radiation field used in the calculation of ^ so that 
Ay = at z = Zb- The material on either side of the 
interface between a clump and the interclump medium 
is in two distinct phases, and we may assume that at 
its outer boundary a clump has a substantial density; 
we take f^(H2) = f^H/2 everywhere in the clump and at 
z — z\^ set ?T.(H2) = 71b- Waves may exist in the interclump 
medium and be partially transmitted into the clump, or, 
as mentioned earlier, may be generated near the inter- 
face by the interaction between the clump and the inter- 
clump medium. Consequently, we assume that the magni- 
tude of the amplitude of the perturbation magnetic field, 
SB = Vbb*, is, at z = Zb, a substantial fraction, /b, of the 
large-scale field, Bq. 

In all models we have taken the mean mass per neu- 
tral particle, ^n, to be 2.3 amu, corresponding to 14% of 
the neutral particles being He and 86% being Hj. As is 
consistent with data given by Savage & Mathis (1979), we 
have assumed 



A. 



1.9 X 1021cm- 



(14) 



where A'h is the column density of Hydrogen nuclei. Also, 
the temperature throughout a clump was taken to be 20 
K. 

For a given setup, an initial boundary value of g{z = 
Zb) = gb, was selected and Eqs. 10, |ll|, and |l2| were 
numerically integrated using an adaptive Gear algorithm 
( pear 1971| ). The value of gh was changed by iteration 
until the inner boundary condition g{z = 0) = was sat- 
isfied. Note that as the total visual extinction (or, equiva- 
lently for a plane-parallel cloud, the column density) for a 
cloud is increased, the clump reaches a maximum size. 
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-^max, and starts shrinking as a larger and larger ther- 
mal pressure (and therefore density) is required to bal- 
ance gravity. Thus, as long as Zb < Zmax, there are two 
values of gh for each setup that satisfy the boundary con- 
ditions. One solution has a relatively flat density profile 
and a small total visual extinction while the other solu- 
tion is more centrally condensed. We deal here exclusively 
with the latter solutions. 

We have considered other models but present full re- 
sults only for models which have a total edge-to-center 
visual extinction of 5 magnitudes since, as discussed in 
R98, it is in the region of Ay of a few that clumps appear 
to begin to contain detected star s, while many dense cores 
may have Ay ^ 5 ( McKee 1989 ) . For our canonical model 
(Model 1) we require a velocity amplitude, V, of 2 km 
sec""'^ and an Alfven speed, va, of 3 km sec~^ at Ay — 2. 
Since the concentration of CO at Ay ^ 2 is very low and 
measurements of GMC clump CO profi les have a width 
of ~ 2 km sec~^ (Williams et al. 1995), observations re- 
quire such velocities to exist well within the cloud. Also, 
we use a wave frequency, of 2 x 10^^^ sec^^; this results 
in relatively strong neutral-ion coupling while keeping the 
wavelength of the perturbing wave less than the size of 
the cloud. In other words. 



Table 2. Summary of Parameters used in Models 



27r 



/ dz/vA 



(15) 



Within the above constraints, we find for our canonical 
model the solution requiring the smallest value of Bo, 
and, since the velocity amplitude of a linear Alfven wave 
is thought to be comparable to but less than the Alfven 
speed, the largest value of /b. However, throughout the 
clump, V < VA, consistent with our assumption of linear 
Alfven waves. 

We present the results of 5 models. Model 1 is for the 
above canonical parameters with R98's depletion case A 
while Model 2 is for case B. Model 3 is the same as Model 1 
but with Lu = Ix 10~^^ to illustrate the effect of a scenario 
with roughly maximum ion-neutral coupling. Model 4 is 
as Model 1 but with an ionization profile given by ^ = 3 x 
IQ-^^Pn"^/^ (see, e.g., McKee 1989 and Myers & Lazarian 
1998). Thus, Model 4 is for a dark region surrounded by 
interclump material. Finally, for completeness, we present 
a model with no turbulence; Model 5 is identical to Model 
1 except with /b = 0. A summary of the parameters of the 
5 models is given in Table |. For all models, rtb = 375 
cm""^ and Bq = 135 fiG. 

5. Results and Conclusions 

In Fig. |l| we present density as a function of visual extinc- 
tion for each of the 5 Models. It is clear that the newer ion- 
ization profiles used in Models 1-3 result in less condensed, 
more extended clumps than the profile used in Model 4. 
In fact, riH is roughly proportional to 1/z in Models 1-3 
while riH goes roughly as 1/z^ in Model 4. In Fig. ^ we 
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Fig. 1. Plot of density (normalized to the density at the 
outer edge of the clump) versus Ay for the 5 Models. The 
solid curve is for Model 1, the dotted curve for Model 2, 
the dashed curve for Model 3, the dot-dashed curve for 
Model 4, and the dash-chain-dot curve for Model 5. 
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Fig. 2. As Fig. |^ but for the perturbing magnetic field, 
5B. 

show plots of SB versus Ay for Models 1-4. Except for 
Model 4, which has hjk^ approaching 0.5 at the clump 
boundary so that Eq. ^ is not satisfactorily satisfied, the 
perturbing field obeys fiux conservation near the surface of 
the clump. In Model 1, in the central region of the clump 
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Fig. 3. As Fig. |i| but for the absolute ion mass fraction, ^. 
The bumps seen on the curves for Models 1-3 are a result 
of the interpolation over Ay. 




Fig. 4. As Fig. |^ but for the magnetic energy flux, F 

k,ULu/\k\^. 
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Fig. 5. As Fig. |l| but for z, the absolute spatial extent of 
the plane-parallel cloud. 



dissipation is rapid enough that SB begins to decrease. 
In order to compensate for the loss of support, the equi- 
librium solution requires a complementary increase in the 
density, as can be seen in Fig. |l|. On the other hand, in 
Model 4, the turbulence is dissipated much nearer to the 
cloud boundary, thus requiring a steeper overall density 
profile. The higher ionization fractions of the case B de- 
pletions result in very little dissipation even in the center 
of the clump for Model 2. N ote that even though observa- 
tions ( Williams et al. 1995 ) suggest that the temperature 
of RMC-type clumps is closer to 10 K rather than the 
20 K used here, thermal support is insignificant except in 
the centre of Model 1, so the effect of a lower temperature 
on the models would be merely to enhance slightly any 
central condensations. 

The ionization profiles used in the Models are shown 
in Fig. ||. The ionization profiles described in Sec. 3 re- 
sult in ^ for Models 1-3 being more than 50 times greater 
near the surface of the clump than in Model 4. However, 
in the center of the clump the ionization fraction drops, 
resulting in more dissipation. Again, this leads to clumps 
that are overall more diffuse but with small condensed 
cores. Clearly, clumps with Ay ^ 5 will have distinct cen- 
tral condensations with n/rib ^ 100 and central fractional 
ionizations of ^ 5 x 10"^. Though dense cores may be 
formed during the fairly rapid collapse (as envisaged by 
Fielder & Mouschovias 1993) of more extended objects 
(i.e. RMC-like clumps) that become unstable, even in our 
equilibrium models we find central cores having densi- 
ties and fractional ionizations similar to those measured 



for dense cores and their envelopes ( 


WiUiams et al. 1995 


WiUiams et al. 199? ; 


Bergin et al. 1999 


)■ 



clumps for Models 1-4. Near the clump center. Model 1 
is nearly thermally supported due to the dissipation of 
the turbulence. The higher ionization fraction for the case 
B depletions used in Model 2 results in less dissipation 
and thus more turbulent support for the clump. Conse- 
quently, as can be seen in Fig. |l|. Model 2 has no central 
condensation and is more extended. Unfortunately, we can 
only speculate about how the depletions of Sulphur, met- 
als, and some oth er species behave in RMC-like clumps 
( Ruffle et al. 1999|) . Thus, cases A and B are merely rep- 
resentative; as can be clearly seen in the figures, the clump 
profiles are very sensitive to the choice of abundances and 
the subsequent fractional ionizations. In addition, com- 
pared to Model 1, the stronger ion- neutral coupling in 
Model 3 results in less dissipation and subsequently the 
clump has little central condensation, as expected. Note 
however that for Model 3 the lower limit in Eq. 15 is not 
adequately satisfied. 

Fig. ^ also shows the effect of external wave genera- 
tion. If the fractional ionization is too low, as in Model 4, 
dissipation occurs close to the surface of the clump. Con- 
versely, if the fractional ionization is too large, as in Model 
2, significant dissipation occurs only at the clump's very 
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centre. Both extremes produce density profiles which lack 
a central condensation. Note that if our externally gener- 
ated wave model is correct, one should not see turbulence 
within a condensed core if there is no turbulence in its 
surrounding envelope. 

In Fig. ^ we present curves which map the visual ex- 
tinction to the spatial extent of the clumps. Clouds with 
larger extents which match the observed 2-3 pc size of 



RMC-type clumps (Williams ct al. 1995) cannot be repro- 
duced within the constraints given in Sect. 4. However, 
observations generally measure the largest linear extent 
of a clump. Thus, since the waves only support the model 
clumps parallel to the large-scale field, it is not surprising 
that the model sizes given here are less than the observed 
sizes. 

Similarly, the models require high boundary densi- 
ties and magnetic field strengths in order for the Alfven 
speed and wave velocity amplitude at visual extinctions 
where CO is abundant to be large enough to be com- 
patible with observed linewidths. For Model 1, rib = 
375 H2 cm^'^. This is rather higher than the typical 
value of n(H2) = 220 cm~'^ given by Williams et al. 
(1995) for RMC-type clumps but, given the uncertain- 
ties, it is within a reasonable range of the Williams et 
al. (1995) value. For Model 1, = 135 /iC, signifi- 
cantly higher than the value of 30 /iG suggested by ob- 
servations ( prlcilcs 1987 ) and expected from robust theo- 
retical arguments (Mouschovias 1987). In order to deter- 
mine whether the values of rib and Bq could be lower and 
still allow model properties to be consistent with observed 
linewidths, we constructed models for clumps with to- 
tal edge-to-center extinctions of 3 magnitudes. The model 
giving V — 2 km sec~^ and wa = 3 km sec~^ at Ay = 2 
had nb = 325 H2 cm~^ = 105 /iG, /b = 0.49, and 
-^max — 0.565 pc; although ?ib and Bq were smaller and 
Zmax larger, the agreement with observations is nonethe- 
less poor. 

Thus, the next step in the modelling of clumps in which 
wave support is important is the inclusion of wave sup- 
port in models analogous to the axisymmetric models of 
magnetically and thermally suported clumps described in 
classic papers by Mouschovias (1976a, 1976b). It is possible 
that the inclusion of magnetic tension, as well as pressure, 
will allow the reduction of Bq to a value more like that 
expected and the construction of models of clumps having 
larger linear extents. 
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